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A STREAMBED-PARTICLE MODEL-STUDY FACILITY 
USING HYDROXYETHYLCELLULOSE 
SOLUTIONS AS A FLUID 


By Neil L. Coleman and Wilbert M. Ellis^ 


ABSTRACT 

Variable-viscosity long-chain polymer solutions will be used to provide 
a wide range of Reynolds numbers in a new facility for experimenting 
with the hydrodynamic forces exerted on streambed particles. Hydrox- 
yethylcellulose solutions tested in the facility proved to have a desirable, 
nearly Newtonian behavior and sufficient viscosity stability. The facility 
can provide the needed range of Reynolds numbers, while maintaining 
flow-velocity profiles typical of streambed conditions. KEYWORDS: hy- 
draulics, hydraulic laboratories, hydraulic modeling, sediment, sediment 
transport, soil erosion. 


INTRODUCTION 

The mechanics of the beginning of streambed- 
particle motion are important in problems of 
stable channel and riprap designs. Many ex- 
perimental studies have been made to investi- 
gate the action of flow-induced forces on a 
streambed particle. In general, where the parti- 
cle is far from the water surface, the action 
of streamflow upon an alluvial particle is a 
Reynolds-similitude phenomenon. For this rea- 
son, experiments have been made using scaled- 
up Reynolds models of a streambed particle. 
In these simulations the character of the local 
flow and the application of forces to the particle 
can be observed in a detail not possible in flume 
experiments. 

In model studies of particles in the silt, sand, 
or gravel size ranges, obtainment of dynamic 
similitude requires operation of the model at 
low-particle Reynolds numbers. Since the model 
particle itself is large, even though the proto- 
type is small, and since the velocities, forces, 
and so forth in the model must be large enough 
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to be measured conveniently, low Reynolds num- 
bers must be obtained by means of working 
fluids of a high viscosity. Furthermore, since 
a wide range of Reynolds numbers is generally 
required in the model, it is advantageous to be 
able to vary selectively the viscosity of the 
working fluid. Early work with a large-scale 
particle Reynolds model indicated that such 
conditions might be accomplished with aqueous 
solutions of a long-chain polymer as a working 
fluid, if a polymer with -sufficiently Newtonian 
viscosity characteristics could be found (5.)^ 

Many long-chain polymers are cellulose com- 
pounds, several of which are used industrially 
as thickeners, sizers, and adhesives. They are 
commercially produced in a variety of types. 
Depending on type, they can be dissolved in 
water to produce non-Newtonian fluids with 
drag-reducing properties or nearly Newtonian 
fluids that have almost stable viscosities over 
a wide range of shear rates. 

Because of their importance in naval engi- 
neering and other applications, the drag-reduc- 
ing polymer solutions have been intensively 
studied for nearly 30 years (d). Perhaps be- 

2 Italic numbers in parentheses refer to items in “Lit- 
erature Cited” at the end of this publication. 
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cause of the emphasis that has been placed on 
drag reduction, polymer solutions with nearly 
Newtonian characteristics have not been used 
as working fluids in model studies where a 
wide range of Reynolds numbers is desired. 

This paper describes a streambed-particle 
Reynolds-model facility designed to operate 
with Newtonian polymer solutions to study lift 
and drag forces on sti-eambed particles, flow 
configuration and pressure distribution around 
particles, and the beginning of particle motion. 
We describe the construction and instrumenta- 
tion of the facility, and the range of experi- 
mental flow conditions attainable by using a 
nearly Newtonian polymer solution. 

THE MODEL FACILITY 

Figure 1 is a schematic of the particle model- 
study facility that has been built, a tunnel that 
recirculates water. It consists of (a) a rectan- 
gular upper conduit within the working section 
of the model, (b) a polymer mixing box, (c) a 
return pipe, and (d) a propeller pump driven 
by a variable-speed device by a 5-hp motor. 
The maximum discharge of the pump is 1,100 
gal/min. 

^ The water tunnel has an auxiliary pump and 
gravity return system (e) for removing fluid 



Figure 1.— Schematic of the model facility: (a) rec- 
tangular conduit, (b) polymer mixing box, 
(c) return pipe, (d) propeller pump, (e) auxiliary 
pump, (f) holding tank, (g) drain, (h) air-vent 
valve, (i) water-filling valve, (j) bleed valve, 
(k) circular opening, (m) packing gland, 
(n) thermometer. 


and storing it in a holding tank (f). This sys- 
tem permits temporary drainage to allow 
access, through appropriate hatches, to the 
working section without wasting the polymer 
solution. For clianging polymer solutions, there 
is ;in external drain (g). Tlie polymer mixing 
box has a removalile lid fitted witii an air-vent 
valve (li) and a water-filling valve (i) with a 
liose fitting. A small lileed valve (j) in the 
upper conduit upstrejim of the working section 
allows withdrawing fluid sainjiles to determine 
viscosity and density. 

To use polymer solutions economically in the 
model facility, an initial solution of a relatively 
high viscosity must be mixed and a complete 
series of exi)oriments performed with this fluid. 
Then, a new fluid with a low viscosity is pro- 


TarLE 1. — Properties of HEC solutions at 77° F 
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1 Last value in each set of measurements was de- 
termined with a falling-ball viscosimeter, other values 
with a rotating-spindle viscosimeter. 

- Each value is the average of 3 determinations made 
at the same shear rate. 
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duced by controlled dilution of the orig*inal 
fluid, followed by another series of experiments. 
The process can be repeated several times, so 
that an extremely wide range of experimental 
viscosities can be obtained with a relatively 
small amount of polymer. 

The upper rectangular conduit (31 ft long, 
0.916 ft wide, and 0.583 ft high) of the model 
facility is constructed of transparent plastic. 
The working section (a) is a bed of Mi-in- 
diameter plastic spheres glued in a close-packed 
configuration to the floor of the conduit, simu- 
lating in scaled-up form a small streambed 
segment. The leading edge of the bed of spheres 
is 22.5 ft from the upstream end of the conduit. 
The bed is 4.6 ft long. At a point on the center- 
line of the bed, and 3.4 ft below the leading 
edge, a 3-in-diameter circular opening (k) in 
the conduit floor allows for the installation of 
equipment for measuring forces on a bed parti- 
cle. A packing gland (m) in the top of the 
conduit, 3.35 ft from the leading edge of the 
bed, allows a total-head tube to be inserted for 
determining the velocity distribution over the 
bed. For monitoring hydrostatic pressure, two 
piezometer taps (not shown) are installed in 
the sidewall of the conduit near the openings 
(k) and (m). A thermometer (n) has been 
installed through the top of the conduit, just 
below the trailing edge of the bed to allow 
fluid temperature to be observed during an 
experiment. 


PROPERTIES OF 
POLYMER SOLUTIONS 

The polymer compound selected for use is 
hydroxyethylcellulose (HEC), which is com- 
mercially available. According to the manu- 
facturer, the compound is nearly Newtonian 
over a wide range of shear rates and relatively 
unsusceptible to viscosity degradation by stir- 
ring or pumping, and it can be completely dis- 
persed in water by stirring before it begins to 
dissolve. Such characteristics insure formula- 
tion of an homogeneous solution, without resid- 
ual clumps of partially dissolved polymer. 

Small samples of HEC solutions at various 
concentrations were analyzed for differences 
in density and viscosity (under varying shear 
rates), to determine the extent to which these 


solutions might exhibit non-Newtonian be- 
havior. These analyses were carried out at a 
standard temperature of 77° F (25° C). The 
results are given in table 1. Each viscosity value 
is the average of three determinations made 
at the same shear rate. 

The solution densities (table 1) were found 
by weighing measured volumes of the different 
solutions. The table shows that the densities 
of HEC solutions at 77° F do not differ ap- 
preciably from the density of water until a 2 % 
concentration is reached. The kinematic viscosi- 
ties in the table were measured with a com- 
mercial rotating-spindle viscosimeter capable of 
exerting eight different shear rates on a sample, 
and also with a commercial falling-ball viscosi- 
meter. In the rotating-spindle viscosimeter, a 
cylindrical spindle with radius r is inserted into 
a small sample container, so that the fluid is 
confined to a restricted annular space with 
radius Ar, small relative to r. The shear rate (??) 
for the instrument can then be approximated 
as a function of spindle speed (w) as 


where Ut is the tangential velocity. 

With the falling-ball viscosimeter, the ball 
falls along the centerline of a tube, so that the 
fluid sample is forced through a small annulus 
at the diameter of the ball. The maximum shear 
rate (v) in the annulus is at the ball surface 
and can be approximated from considerations 
of the velocity distribution along any radius of 
the annulus, and shown as a function of the 
ball Reynolds number. 



where Rt, is the ball radius, W the ball velocity, 
and V the kinematic viscosity of the fluid. 
Apparently (table 1), the falling-ball viscosi- 
meter exerts substantially higher shear rates 
than does the rotating-spindle viscosimeter. 

The non-Newtonian behavior of the HEC 
solutions is not pronounced: kinematic viscosity 
for the six solutions in table 1 decreased about 
34^0 over an average 100% increase in shear 
rate. The variation of kinematic viscosity with 
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shear rate for these polymer solutions can be 
represented by 



where vi is the kinematic viscosity at the unit 
shear rate, rji. For a perfectly Newtonian fluid, 
the exiDonent N would have the value zero. For 
the polymer solutions, N varies with concentra- 
tion, as shown in figure 2. The relative-viscosity 
function expressed by equation 3 is depicted in 
figure 3. 

Evidently, there are two ways of representing 
the apparent viscosity of a polymer solution 
used as a working fluid in the model facility. 
The first is to use the unit shear rate (?;,) for 
a particular solution, at a particular fluid tem- 
perature. However, this is not realistic with 
regard to actual shear rates expected in the 
model. The second method is to use the nearly 
constant (quasi-Newtonian) viscosity exhibited 
by the polymer solutions at high shear rates. 
We adopted this method to study the flow 
conditions in the model facility, since the high 
shear-rate viscosity of a given working fluid 
can be easily obtained with the falling-ball 
viscosimeter. 

PERFORMANCE OF THE 
FACILITY 

We tested the series of HEC solutions for 
which temperature-viscosity functions are given 
in figure 4, to verify that the polymer solutions 



Figure 2.-Variation of the exponent N with polymer 
solution concentration. 


would maintain constant, viscosity duri 
pumping and to define the velocity distributio 
and the range of experimental conditions th 
could be generated in the facility. 

The viscous stability of the HEC solutio: 
was verified by comparing the temperatui 
viscosity function immediately after mixii 
and after several hours of experimental pum 
ing. Figure 5 shows the comparison for a 0 3 ' 
HEC solution. 

Because of the wide range of viscosity of tl 
HEC solutions, the flows generated in tl 
facility can be expected to display time-mea: 
velocity distributions ranging from the par 
bolic distribution of a completely viscous flo 
to the logarithmic profile of a turbulent, rou| 
conduit flow. These velocity distributions eg 
all be accommodated (7) on a single grap 
plotted in the coordinates 



where U is the time-mean velocity, y the di; 
tance from the velocity measurement to tt 
virtual origin of the profile, [/,„ the shear v( 
locity at the virtual origin, « the Karman coi 
stant, and S the boundary-layer thickness, take 
as the distance from the point where the max: 
mum velocity [/,„ occurs in the conduit to th 
virtual origin. The ratio aU/U-^: is the velocit 



Figure 3.—.Relative-viscosity function for HEC solu 
tions. 
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Figure 4. Viscosities of working fluids for solutions of various nominal concentrations. 


shift caused by boundary roughness, and the 
function (H/k) [(o(4//8) ] is Coles’ wake func- 
tion iU) , in which n is the coefficient of wake 
strength. The function fiUty/v) is the “stand- 
ard” velocity-profile function that can be de- 
rived from simple dimensional analysis or from 
mixing-length theory (5). 

Velocity profiles were measured in the model 
facility with an appropriately designed and 
calibrated total-head tube (2) and a differ- 
ential pressure transducer. The total-head tube 
was inserted through the packing gland (fig. 1, 
m) . It was constructed so that its tip would 
traverse a vertical located 3.2 ft from the lead- 
ing edge of the test bed, or just ahead of the 
transducer insertion port (fig. 1, k). In making 
the velocity measurements, we took the virtual 
origin of the profile as one-fifth of a sphere 
diameter (0.2Z) or D/5) below the top plane 
of the bed, as recommended by Einstein and 
El-Samni (5). 

The velocity profiles were plotted in U/Um 
(UmV/v) coordinates, and U* and £i.U/U» were 


estimated, using an adaptation of the cross-plot 
method (8). 

In figure 6, the velocity shift Af//f/* is plotted 



Figure 5. — Viscosity-temperature curve for a 0.3% 
polymer solution (a) after initial mixing, (b) after 
7 hours of pumping during experiments. 
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Figure 6. — Velocity-shift function, 
against the bed-shear Reynolds number: 


AU u..D\ 

— =5.6 log 4.2, (6) 

which is identical to that for sand roughness, 
except the constant is 4.2 rather than 3.6, indi- 
cating that a bed of closely packed spheres is 
slightly smoother than a bed of natural sand 
particles. The above results agree generally 
with those of other experimenters (10). 

The velocity profiles are plotted in figure 7. 
The functions included in the figure are, except 
for the introduction of aE//Z 7«, the usual 



for the viscous sublayer, and the Karman- 
Prandtl law of the wall, 


R,s= 


V,I> 


(5) 


U 

Tl. 


2.303 


log 




+4.9, 


(8) 


Figure 6 shows that bed roughness begins to 
affect the flow at C/..J9 / v equal to approximately 
4. The transition regime ends, and the fully 
rough flow begins at about 50. In the fully 
rough regime, the function is 


for the inertial region. For plotting figure 7, 
the Karman constant, k, in equation 8 was 
taken as 0.41. The dash curve in the figure is 
the usual empirical transition function between 
the viscous and inertial regions. The velocity 



V 


Figure 7. — Typical velocity profiles. 
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profiles labeled A in the figure are characteris- 
tic of totally viscous flow and have the equation 



which is expanded to 



for plotting in figure 7. These profiles are 
asymptotic to the viscous-sublayer function. 
The profiles labeled B are typical of flows 
where the displacement Reynolds number, 



is less than 72r) (,S’). 'riiese flows have a strong 
wake component and no logaritlimic region. 
Profiles lal)eled C are tyi)ical of flows with 
Rgsi greater than 725, with significant intensity 
of turbulence in the inner part of the conduit. 
These profiles sliow botli a logarithmic part and 
a wake region. As conduit tui'l)ulence increases, 
the wake component decreases (-0. Jincl pi'o- 
files like D appear. These profiles have a loga- 
rithmic region l)ut no wake region. 


The bed-shear Reynolds number (R*) is the 
independent variable on which the velocity re- 
duction in figure 6 depends. This Reynolds 
number is an indicator as to whether the flow 
over the bed is hydraulically smooth or rough, 
and whether or not a particle on the bed is 
experiencing viscous or turbulent approach 
flow. Thus, R* is one type of particle Reynolds 
number, a second is the conventional form 


where U» is the approach-flow-velocity vector 
acting along a line through the particle center. 
This Reynolds number characterizes the local 
flow around an isolated streambed particle and 
serves as the primary independent similitude 
parameter in lift and drag coefficient studies. 

The joint variation of R.i, and Rc is plotted in 
figure 8. Also included for comparison are 
the ranges of R,. obtained by Coleman (3) and 
Aksoy (i), and the ranges of R* obtained by 
Watters and Rao {9). Apparently, the facility 
described here has broader experimental po- 
tential than those used in the other model 
studies. The capability of operating over a 
range from totally viscous flows to fully tur- 
bulent rough boundary flows is particularly 
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advantageous for streambed-particle model 
studies. 


CONCLUSIONS 

By using^ a water-tunnel facility and HEC 
solutions with a wide range of viscosities as 
working fluids, researchers can conduct stream- 
bed-particle model studies over a wider range 
of Re5m.olds numbers than obtained previously. 
Particle Reynolds numbers from 1.0 to 10,000 
can easily be obtained, and flow conditions can 
range from completely viscous to completely 
hydraulically rough. 
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